A quantum dot implementation of the quantum NAND algorithm 
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We propose a physical implementation of the quantum NAND tree evaluation algorithm. Our 
approach, based on continuous time quantum walks, uses the wave interference of a single electron in 
a heirarchical set of tunnel coupled quantum dots. We find that the query complexity of the NAND 
tree evaluation does not suffer strongly from disorder and dephasing, nor is it directly limited 
by temperature or restricted dimensionality for 2-d structures. Finally, we suggest a potential 
application of this algorithm to the efficient determination of high-order correlation functions of 
complex quantum systems. 

PACS numbers: 03.67.-a,73.23.-b 



Recently a new quantum algorithm was discovered that 
evaluates a binary tree of incommeasurate conditions via 
the logical NAND gate in a time faster than the prov- 
ably best classical algorithm 0, Q, H, 0, @|. One con- 
struction of this algorithm relies upon a continuous time 
quantum walk between coupled bound states, where the 
"answer" to the posed query is determined by the re- 
flection or transmission of a single-particle wave-packet 
through an idealized, complex scattering device 0, @]. 
Mesoscopic transport in semiconductor heterostructures 
have ideal properties for creating a device to implement a 
continuous time quantum walk. Their long phase coher- 
ence times and lengths for electronic transport allows for 
novel studies of the electronic properties of small-scale 
structures @, Modern fabrication techniques allow 
for complex structures with robust controls and fast ma- 
nipulation [^, [T3, 11, EH) an d the key features of the 
continuous time quantum walk approach (a Fano-like in- 
tereference effect with a non-trivial set of coupled bound 
states) have been experimentally demonstrated (l3j . 
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In this paper we propose an approach to implementing 
NAND trees by mimicking the logical structure of the 
Hamiltonian formulation with a physical arrangement of 
quantum dots. Using a recursive Green's function ap- 
proach inspired by Ref. 0, we find the quantum speed-up 
of the quantum NAND algorithm — that it takes a time 
0{\fN) for N input bits — is maintained if the phase de- 
coherence rate l/r^ @] and the disorder (detunings and 
tunnel couplings) be lower than t/y/N, where t is the 
characteristic strength of interdot tunnel coupling t and 
N is the number of bits in the input to the NAND tree. 
One key feature of the present work is the use of a large 
physical space (O(N) quantum dots are used) to allow 
for efficient coupling to the complex quantum system (the 
"oracle" ) to which the query is posed. We develop appro- 
priate techniques such that temperature and restricted 
dimensionality lead to at most a O(logiV) overhead in 
the evaluation of the algorithm. We conclude with a 
discussion of the "quantum-ness" of the algorithm and 
consider how such a device can be used to measure many- 
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FIG. 1: The logical NAND tree with N = 8 inputs in the 
top of the tree. Nodes are labeled from 1 to 15. Each node 
can be realized as a quantum dot, with the links between 
nodes represented by tunnel coupling between linked dots. 
Detunings bi applied to dots i + N at the "top" of the tree 
set the local potential to be off of the Coulomb blockade peak 
for dot i + N ii bi 0. The conductance of the dot i — at 
the bottom of the tree, near zero bias (E = 0), encodes the 
result of the NAND operation. 



body properties of complex quantum systems. 

Quantum dot implementation: The NAND tree quan- 
tum walk is described by the transmission coefficient of 
a series of coupled states, shown in Fig. [T] Each "state" 
|i) in the graph corresponds to the lowest orbital state 
of a quantum dot in the Coulomb blockade regime, with 
i = l labeling the bottom quantum dot at level k = 
and i = j + 2 k — 1 labeling the jth dot at the A:th level. 
We assume a large orbital level spacing a ^> t, where t 
is the average tunnel coupling in the tree of dots. The 
coupled dot Hamiltonian can be written: 



H = 2^ ei\i)(i\ - t 2l \i) (2t| - t 2l+l \i) (2i 



ll+H.c. (1) 



where \i) = c\ \rii) describes the addition of a (spin-less) 
electron on the ith dot near the charge transition from 
Tii — ► Tii + 1. Each dot can have a slight detuning from 
the Fermi energy of ej , and each dot i is tunnel coupled 
to dots 2i and 2i + 1 with a tunnel coupling tn, t-a+i ~ t, 
forming a physical binary tree structure. 

At the top of the tree is an oracle (O), a quantum 
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system with electron number operators hn coupled to 

the binary tree via V = AX^o^ - )^'^ + ^X* + ^1- 
Unlike previous work, we focus on a coupling that may 
be easy to implement physically: for example, the oracle 
could be a mesoscopic device capacitively coupled to the 
top of the tree. We will proceed with the case of an oracle 
in an eigenstate of the {hio} operators with eigenvalues 
{hi} each taking values or 1. This sets the detuning of 
dots at the top of the tree to ti+N = {— YhA. 

The algorithm progresses by the injection of an elec- 
tron wavepacket into the bottom of the tree. This probes 
the Green's function of the bottom quantum dot, which 
depends recursively on the Green's functions for the left 
and right subtrees. For wavepackets with small energy 
spreads AE < t/y/N, the successive, coherent scattering 
events up and down the tree lead to either total trans- 
mission or total reflection of the wavepacket, where the 
result (transmission or reflection) is the final "answer" to 
the NAND query for an input specified by the bits 6, . A 
voltage probe in the linear response regime provides such 
an energy probe for the system: the conductance for elec- 
trons passing through the bottom of the tree encodes the 
result of the tree. To show that our approach works in 
the presence of dephasing and disorder, we develop a re- 
cursive description of the scattering Green's function in 
the quantum dot tree in the spirit of previous work Q • 

Green's function approach: We now describe an ap- 
proach to understanding the NAND tree evaluation using 
a recursively defined Green's function for single-particle 
propagation through the tree. If we consider a "tree" 
segment with a single state, 1, at the bottom, we can 
define the projectors Pi,Qi,Q r , where Qi projects into 
the subspace associated with the left side of the tree (all 
states connected to dot 2) and Q r similarly projects into 
the subspace associated with the right side. We can write 
the inverse Green's function exactly [3]: 



G^(E) 



E + ij- P 1 HP 1 



3=1 ir 



E + i'y- QjHQj 



QjHPi 



= E + i 1 -e 1 -\t l \ 2 Gi{E)~\t r \ 2 G r {E) (2) 

This gives a recursive relation for the projected Green's 
function Gi (E) for a state i in terms of the Green's func- 
tions for the subtrees connected to states I = 2i and 
r = 2i + 1. The term 7 = l/r^ is a measure of the 
charge-based dephasing the electron may be expected to 
encounter. 

At the top (largest fanout) of the tree, we can write the 
single-state Green's function G t ° v (E) = E+ l . which 
truncates the recursion relation. We remark that this 
recursion is only successful because the tree contains no 
loops; otherwise, Eq. ([2]) would fail. In the disorder-free 



case, ti = t for all i and e, = except for dots at the top 
of the tree with input bits bk ^ 0, as described above. 

Setting energy to units of t, we consider the Green's 
function at the top of the tree for the three different in- 
puts. We set Si to be the sign of the detuning to illumi- 
nate the role the alternating sign (— ) l plays. We assume 
7, E -C 1 < A. For dot j, tunnel coupled to 2j and 2j + l, 
we get (to 0(t/A,j/t)): 



Gj(E) 



-E/2-ij/2 

—E — i"f 
1 



hj = 0, &2j+i = 
hj = 0, &2j+i = 1 
hj = 1,&2j+i = 1 



(3) 



As the Si alternate such that each paired "1" input has 
S2j = — S2j+i, the third term has no energy shift. We 
now identify Green's functions of the forms 



G- 



1 



aE + «7/3 



Gi» = -(aE + ij/3) (4) 



where the "0" and "1" indicate the expected NAND re- 
sults for the inputs. Any dependence on A is removed 
by the alternation of Sj for paired "1" inputs. 

We can include effect of disorder in the system, both in 
variations in ti and in variations in e^. Both are likely to 
occur due differences in fabrication, low frequency noise 
on electrostatic gates (such as those used to tune the 
system to the multi-dot regime) and charge-fluctuators 
in the system. For the three different inputs (00,01,11) 
at the top of the tree: 



0,0 



-E 



jei - i-/) 2 tj + (e r - ^ftl 
tf + ei4 - i 7 (t? + tl)f 



1,0 

1,1 



-E 



ei-tf e r -ig 

_L + e ' r ~ ' 7 



E+(i-y+(sit 2 + M?)/A - e,-) 



(5) 

(6) 
(7) 



These have the same form as the disorder free results of 
Ref. so long as ti,t r > t, where t is now the minimum 
tunnel coupling. More generally, we will assume U and 
ti to be drawn from a normal distribution with means 
t, and standard deviations er t , o~ e . We can remove any 
explicit A dependence as long as at/ A <C 1. 

We now show that two subtree's described by Green's 
functions of the above type lead to the appropriate type 
Green's function at the joining of the two subtrees, re- 
cursively implementing the NAND tree. Consider the 
general map for two subtrees' Green's functions Gn-^j) 



and G 



r(=2j+l), 



parameterized by ai r ,(3i r . Gj for the 



cases 00, 10, and 11 are given to O^/t, £j/t, t/A): 



3 



1 



1 



aiE + i-yfli a r E + ijf3 r 
-(aiE + i"fPi), -{a r E + ij0 r ) 



otitic 



(t 2 r ft + tf0 r 
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t 2 r Pl + tjP r 



(1 + tfai + t 2 r a r )E + + iffy + t 2 [3 r + iej/j) 



(8) 
(9) 
(10) 



The overall map implements a NAND gate by identifica- 
tion of "0" and "1" via Eq. H 

For a tree consisting entirely of the "most dangerous 
case" (an input of the form 1011) subtrees of Ref. |2j, we 
find a — > 2 fe / 2 , and f3 — > 2 k / 2 + i J2j e j/7, where the sum 
is over the set of all dots where both subtrees evaluate 
to 1. The first term indicates that the Green's func- 
tion derived is only appropriate within an energy range 
t/2 n / 2 = t/VN about E = 0. This leads to the condi- 
tion that the electron wavepacket be spread over a time 
~ /t, and sets y/W scaling of the evaluation time. For 
random, uncorrelated noise, the second term (with O(N) 
terms) has a root-mean-square of cr £ y/~N, and leads to an 
energy shift by that amount from the desired resonance 
at E = as well as a dephasing r )2 n l 2 . This indicates 
that 7,cr e <C t/y/W form the requirements for operation 
of the device. Less difficult is the bounding requirement 
U > t. We compared these results with numerical simu- 
lations and find consistent behavior for tree sizes up to 
N = 128. A sample with N = 32 is shown in Fig. O 

Measuring with a quantum dot: We now detail how a 
quantum dot coupled to a linear voltage probe can mea- 
sure the result of the NAND tree well below the tem- 
perature limit of the leads. This is a generalization of 
sub-temperature spectroscopy used in double quantum 
dot transport measurements [9(. 

We write the Green's function, including the self en- 
ergies £f r for left and right leads coupled to the dot via 

tl,r: &S 



Go\E) =E~e Q 



M 2 Gi(£) (11) 



By working well above the band-bottom of the leads, we 
can write T = Ti + T r with T li7 . = 2t 2 l r v = -2iE^ r , 
where v is the density of states near the Fermi energy. 
The Meir formalism [l5( for the transmission Ti r {E) = 
Tk[T,Gi ) or r Gf ] yields 



f lr (E) 



r 2 /4 



(t 2 Im[Gi(£)] - r/2) 2 + (E- t 2 Be[Gi(E)} - e ) 2 

(12) 

The conductance follows from the Landauer formula, 
Gw = e 2 /h J f(E - Ef)fi r (E)dE, where f(E) is the 
Fermi distribution. 

A few remarks are in order. First, a pole at E = 
for Gi(E) leads to zero conductance, while a "l"-like 




FIG. 2: Simulations for N = 32 with outputs "0" 
(blue, purple) and "1" (red, orange). The dependence of con- 
ductance on eo (left) and E (right) are shown. Disorder 
parameters 7 = <r e = a t — 0.03t (red, blue) and O.lt (or- 
ange, purple) are used. 



Green's function (proportional to energy), leads to a peak 
at eo = 0, thus letting the existence or absence of trans- 
port indicate a "1" or "0" result. Second, when dephas- 
ing (7) is weak, the transmission function is a Lorentzian 
with a width at most r/2. When convolved with the 
derivative of the Fermi function at k^T > T, the edges 
of the Fermi function are effectively damped, leading to 
sub-temperature energy resolution (limited only by T). 
To see these effects, we simulate the expected conduc- 
tance as a function of eo and E at T — in Fig. [2] for an 
input of N — 32 bits. Cases resulting in "0" and "1" are 
both shown. For this scenario, we set T/t = 0.1 < 2 -5 / 2 
and A/t = 10 > 2 5 / 2 . The temperature dependence 
can be mitigated by further reductions of T. However, 
disorder, as described below, may limit the achievable 
resolution. 

Two-dimensional layout: One difficulty encountered 
for making large NAND trees using this approach is the 
increasing perimeter of the tree as a function of tree depth 
(N versus log N) . For sufficiently large trees, the spacing 
between successive levels k and fc+1 of the tree must grow 
as 2 fc /( d_1 ) for a d-dimensional layout. Furthermore, tun- 
nel coupling exponentially decreases as the spacing be- 
tween quantum dots increases, potentially limiting this 
approach to small trees. 

As an alternative, "dummy" levels of the tree can be 
inserted: these are states coupled only to a single dot 
going up or down the tree (Fig. [3^)- Each acts as an 
inverter, where the "missing" state is equivalent to the 
absence of a resonance in one "leg" of the tree, which 
is the "l"-type Green's function. Specifically, a state of 
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FIG. 3: (a) An inline quantum dot act as an "inverter" — 
the non-existant right inputs to each dot act as a virtual "1" . 
Two inline dots move the bit in the top input to the bottom 
output, (b) An array of tunnel coupled quantum dots in an H- 
fractal layout, with pairs of inverters connecting more distant 
quantum dots. The number of inverters between dots doubles 
every two levels of the tree. 



either "0" or "1" type parameterized by (a, (3) becomes, 
after coupling through two inverters, (1+a, 1+/3). Going 
a distance d (through 2d inverters) gives 



(d + a, d + (3) 



(13) 



As an explicit construction, we consider the H-fractal 
(Fig.[3jD) for the layout of quantum dots. The H-fractal is 
a binary tree structure comprised of an "H" with the two 
sides the same length as the middle bar. [l6[ This has 
the benefit that the distance between successive levels of 
the binary tree doubles for every two levels. A "worst 
case" trees (again, with inputs 1011) has the mapping 



K+2, Pk+2) = (2 fe/2 + 2a, 2 fc / 2 + 2/3) 



(14) 



Where the distance goes as 2 fc / 2 at the fcth level (with 
n the bottom of the tree). A tree comprised entirely of 
worst-case trees has 

a n = (3 n = 2 n ' 2 + 2{2 n ' 2 - 1 + 2(. . . (1) . . .)) (15) 
< n2 n/2 = log 2 (N)VN (16) 

This indicates that a 2D architecture leads to an over- 
head of log(AT) for the evaluation of the tree. This still 
performs better than the best classical algorithm. As a 
final remark, we can selectively remove one of the pair of 
dots to produce a "NOT" operation inline with the fol- 
lowing NAND operation by opening the tunnel coupling 
between the pair of dots. This allows the tree evalutor 
to evaluate arbitrary boolean functions, not just NAND 
trees. 

Practical limits: Taking values 7 = 0.1/ieV, t = 
100/^eV, a = 1 meV, T — O.ltteV, a e = a t = liteV, and 
fc&T = 2/ieV (20 mK dilution refrigerator electron tem- 
perature), we find the limiting algorithm size N ~ 2 13 , 
where the limit is due to the disorder in the system. Us- 
ing the H-fractal layout, this corresponds to a mesoscopic 



region of area a x 3 13 < 0.02 mm 2 for dot center-to- 
dot center spacing a ~ 100 nm. This includes appro- 
priate space for electrostatic gates or etched barriers. 
Phase coherent transport over such large regions have 
been demonstrated in a number of experiments Q. The 
expected evaluation time is of order 10/r < 60 ns (to 
generate a differential signal of ten electrons, sufficient 
for fast detection). 

We can use this NAND tree evaluator as a subsystem 
in a larger processing unit. The best known classical al- 
gorithm evaluates in jV 0,753 time. Using the quantum 
device as a sub-tree evaluator of a quantum oracle, a 
problem of size 2 13+fc can be evaluated in 2 6 - 5+a753fc < 
29.79+0.753'° tml6j as su btrees of size < 2 13 can be cal- 
culated sequentially using the quantum subsystem as a 
hardware accelerator. Better hybrid quantum-classical 
approaches may exist, but require further investigation. 

Conclusion: This paper illustrates a fundamental and 
heretofore unrecognized quality of the quantum NAND 
algorithm. In particular, the entire NAND quantum walk 
can be done with a single-electron wavepacket and ap- 
propriate resonances (quantum dots). Other approaches 
using classical waves, such as light in coupled cavities, 
would in principle have the same Green's function. Thus, 
from a quantum information perspective, the interesting 
part of the algorithm lies in the oracle: the quantum me- 
chanical system which "sets" the inputs. We heretofore 
considered the oracle to be in an eigenstate of the num- 
ber operators nto- If we describe the boolean function 
evaluation of the tree by /(&o, ■ ■ ■ , frjv-r) = or 1, an ar- 
bitrary oracle state leads to an expected result of the tree 
of (f(noo, . . . , fiN-io))o- For example, a two-level tree 
evaluates (n e>«re> + «20«-30 - ^oo^io«20«3C>), a non- 
trival, high-order correlation function of the oracle. We 
now have the outstanding task of identifying interesting 
quantum systems to use as an oracle. 
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